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ABSTRACT 

X-ray ionized reflection occurs when a surface is irradiated with X-rays so intense that its 
ionization state is determined by the ionization parameter £, (x F/n, where F is the incident 
flux and n the gas density. It occurs in accretion, onto compact objects including black holes 
in both active galaxies and stellar-mass binaries, and possibly in gamma-ray bursts. Compu- 
tation of model reflection spectra is often time-consuming. Here we present the results from 
a comprehensive grid of models computed with our code, which has now been extended to 
include what we consider to be all energetically-important ionization states and transitions. 
This grid is being made available as an ionized-reflection model, REFLION, for XSPEC. 
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1 INTRODUCTION 

The process of backscattering and fluorescence of X-rays from 
cosmic matter known as X-ray reflection is commonly observed 
throughout astronomy. We are concerned here with the situation 
where the X-ray irradiation is so intense that it determines the ion- 
ization state of the material. This occurs in accreting X-ray sources 
such as Galactic Black Hole Candidates (BHC) and Active Galac- 
tic Nuclei (AGN), and may also be relevant to Gamma-Ray Bursts 
(GRB). 

The irradiated gas is Compton thick, which requires that 
Comptonization of the radiation as well as the necessary atomic 
physics be dealt with. Results of such computations have been pub- 
lished over the last ten years, first for constant-density atmospheres 
(Ross & Fabian 1993; Zycki et al. 1994; Ross, Fabian & Young 
1999) and more recently for density structures in hydrostatic equi- 
libium (Nayakshin, Kazanas & Kallman 2000; Ballantyne, Ross & 
Fabian 2001; Rozanska et al. 2002; Mauche et al. 2004). Reason- 
able agreement is found between most of these results (Pequignot 
et al. 2002). 

Most methods are computationally very time consuming, so 
only limited regions of parameter space have been explored. Our 
method runs relatively rapidly, on the other hand, when the illu- 
minated gas is assumed to have constant density. Such models are 
still useful for several reasons. The overall geometry is still not 
understood for many sources. Even when a simple accretion-disc 
structure is assumed, detailed calculations of hydrostatic equilib- 
rium are subject to uncertainties in the total thickness of the disc 
(and hence the vertical component of the local gravitational field), 
the boundary condition (either pressure or height) at the base of 
the illuminated layer, the external pressure exerted by the corona. 
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the geometry of the illumination, the influence of magnetic fields, 
etc. Furthermore, Ballantyne et al. (2001) found that models for re- 
flection by atmospheres in hydrostatic equilibrium could be fit by 
diluted versions of constant-density reflection models. 

We have now assembled an extended grid of models for which 
the ionization parameter, irradiating spectral index and iron abun- 
dance are free parameters. At the same time the original code has 
been expanded to include further ionization species. All abundant 
species and their important transitions are now included. The pur- 
pose of this paper is to present example spectra from across this 
grid which can be used as a guide to the interpretation of AGN, 
BHC and GRB X-ray spectra. 



2 METHOD 

2.1 Radiative Transfer 

Our calculations extend and improve on the models discussed by 
Ross, Fabian & Young (1999). Now the illuminating radiation is 
assumed to have a cutoff power-law spectrum, 

Fe = AE-''+^ exp(-£/Sc), (1) 

that extends to higher photon energies. Different values of the 
photon index P are treated, while the cutoff energy is fixed at 
Ec — 300 keV. The amplitude A is chosen so that the total il- 
luminating flux, Ftot ~ J Fe dE, corresponds to a desired value 
of the ionization parameter. 



47rFtc 



(2) 



The illumination is incident on the upper surface of a plane- 
parallel slab representing the top layer of an optically thick medium 
such as an accretion disc. The slab has hydrogen number density 
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riH = 10^^ cm~^ and Thomson depth tt ^ 5. (Larger Thom- 
son depths are treated for larger values of ^, since the illumination 
can penetrate to greater depths when the gas is more highly ion- 
ized.) No net flux is allowed to enter the slab from below, so that 
the emergent radiation is due entirely to reprocessed illumination 
("reflection"). Additional radiation from beneath the surface layer 
could be approximated by adding a soft blackbody spectrum to the 
calculated reflection spectrum (see Section 4). 

As described by Ross & Fabian (1993), the penetration of the 
illuminating radiation is treated separately from the transfer of the 
diffuse radiation produced via scattering or emission within the gas. 
In order to extend the treatment of the diffuse radiation to photon 
energies above 100 keV, Compton scattering is treated using the 
Fokker-Planck equation of Cooper (1971): 

d \ ,^ „^^A f , ^ 9w\ 
dE) 
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(3) 



where n is the photon occupation number, E is the photon en- 
ergy. Tic is the free-electron number density, and 8 — kT. Here the 
Fokker-Planck equation of Kompaneets (1957) has been modified 
by the term 



a{E,e) = 



l+/3{e)/{l + 0.02E) 
1 + 9 X lO-'^E + 4.2 X 10-<^E^ 



(4) 



(5) 



with E in keV and with 

Equation (3) is applicable for E < 1 MeV. The diffuse radiation 
field is taken to satisfy the steady-state Fokker-Planck/diffusion 
equation. 
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where z is the vertical height within the slab, Je is the spectral 
emissivity, k — ka + i^kn is the total opacity (per volume), ka 
is the absorption opacity, and kkn is the Klein-Nishina opacity 
for Compton scattering. We apply equation (6) to radiation with 
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< £ < 10 keV. The energy resolution is E/AE 



50 



throughout the 0.1-10 keV spectral range (except for a resolution 
of 70 around the iron Ka lines), with somewhat lower resolution 
outside that range. 



2.2 Atomic Data 

As the radiation field is relaxed to a steady state, the local temper- 
ature and fractional ionization of the gas are found by solving the 
equations of thermal and ionization equilibrium, as described by 
Ross & Fabian (1993). For this non-LTE calculation, we have ex- 
tended the number of elements and the range of ions treated, and 
we have updated much of the atomic data employed. In addition to 
fully-ionized species, the following ions are included in the calcu- 
lations: C III-VI, N III-VII, O III-VIII, Ne III-X, Mg III-XII, Si IV- 
XIV, S IV-XVI, and Fe VI-XXVI. 

Photoionization cross sections for all subshells of the ions 
treated are calculated from the fits of Verner & Yakovlev (1995). 
Rates for direct collisional ionization and excitation autoionization 
are taken from Amaud & Rothenflug (1985) and Arnaud & Ray- 
mond (1992). For carbon through sulfur. Auger ionization is treated 
as described by Weisheit & Dalgamo (1972). For iron, we employ 
the probabilities for multiple Auger ionization calculated by Kaas- 
tra & Mewe (1993), and a generalization of the method of Weisheit 



(1974) is used to calculate the resulting effect on the ionization bal- 
ance. 

For all elements except iron, total recombination rates as func- 
tions of density and temperature are derived from the tables of 
Summers (1974). These rates include two effects that are impor- 
tant at high densities: the increase in the recombination rate at low 
temperatures due to three-body recombination and the reduction 
in the dielectronic recombination rate at high temperatures due to 
collisional ionization of the highly excited states that follow radia- 
tionless recombination (Burgess & Summers 1969). At low densi- 
ties, these rates are in reasonable agreement with more modem fits 
for combined radiative and dielectronic recombination (Arnaud & 
Rothenflug 1985). For iron, rates for radiative and dielectronic re- 
combination are calculated using the fits of Amaud & Rothenflug 
(1985) and Vemer & Ferland (1996). Three-body recombination 
of iron is neglected, since it is unimportant at the densities under 
consideration (Jacobs et al. 1977). Rates for radiative recombina- 
tion directly to ground levels are calculated using the Milne relation 
(e.g., Bates & Dalgamo 1962). 

The most important emission lines for each ion are treated 
in the calculations. Since the illuminated atmosphere is thick and 
dense, resonance lines are assumed to be optically thick ("case B"), 
greatly reducing the number of important lines. A total of 274 mul- 
tiplets are included; individual lines within a multiplet are not re- 
solved. Line energies and oscillator strengths for resonance lines 
are taken from the results of the Opacity Project (Vemer, Verner 
& Ferland 1996). Excitation and line energies of nonresonance 
lines, as well as collision strengths for all lines, are derived from 
the CHIANTI database (Dere et al. 2001). Parameter values miss- 
ing from these databases are taken from Kato (1976), Raymond 
& Smith (1977), Gaetz & Salpeter (1983), or Landini & Mon- 
signori Fossi (1990), as necessary. The iron Ka lines treated are 
the recombination lines of Fe XXVI and Fe XXV (near 7.0 and 6.7 
keV, respectively) and the fluorescence lines of Fe VI-XVI (near 
6.4 keV). Kq fluorescence of Fe XVII-XXII is assumed to be sup- 
pressed by autoionization (Auger effect) during resonance trapping 
(Ross, Fabian & Brandt 1996). Recombinations to excited levels 
of Fe XVII-XXII are assumed to result ultimately in 3s 2p line 
emission (Liedahl et al. 1990). 

Resonance lines require special attention. Resonance-line 
photons can avoid destruction during resonance trapping by two 
mechanisms: Compton scattering out of the narrow line core be- 
tween resonance scatterings and reemission in the far line wings 
during resonance scattering. The fraction avoiding destmction is 
taken to be 



/ = 



(KT/«:)-Pcont + Pose 
+ -Pcoll + -Pcsc 



(7) 



Here kt is the Thomson opacity, k is the total continuum opac- 
ity at the line energy, Pcont is the probability per resonance scat- 
tering that the photon undergoes a continuum process (absorption 
or Compton scattering), and Pcoii is the probability per resonance 
scattering that collisional deexcitation occurs instead of reemission 
(Hummer 1968). The escape probability Pcsc due to reemission in 
the far line wings (Hummer & Rybicki 1971) is calculated using the 
approximation for the K2 function given by HoUenbach & McKee 
(1979). A Doppler line profile is assumed for all resonance lines 
except the Fe Kq line, where a Lorentz profile is assumed due to 
the importance of the damping wings. Line photons that avoid de- 
struction during resonance trapping are added to the continuum by 
means of the local emissivity. 

Elements lighter than iron are assumed to have solar abun- 
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Figure 1. Structure as a function of Thomson depth for the isolated slab 
illuminated on one side. The upper panel shows the temperature, while the 
middle and bottom panels show the ion fractions of oxygen and silicon. 



dances (Morrison & McCammon 1983), while the abundance of 
iron is set to a factor ^fc times its solar value. Models are gener- 
ated with underabundances as well as overabundances of iron. On 
the one hand the X-ray spectra of some Seyfert galaxies appear to 
require high abundance of iron, for example MCG-6-30-15 (Fabian 
et al. 2002) and IH 0707-49 (Boiler et al. 2002), and Shemmer et 
al. (2004) have reported a strong correlation between metallicity 
and accretion rate in a sample of luminous quasars. On the other 
hand, Nomoto et al. (1997) have found that Type II supemovae 
produce the following average abundances relative to solar values: 
Ac = 0.2, An = 0.01, Ao = 1.3, ^Ne = 0.9, Ams = 1-2, 
Asi = 1.1, As — 0.7, and Afc ~ 0.35. In that case, oxygen, 
neon, magnesium and silicon are near solar abundances, while the 
abundance of iron is markedly reduced. 



2.3 Test Case: an Isolated Slab 

Recently Dumont et al. (2003) have criticized calculations, such as 
ours and those of Nayakshin et al. (2000), that employ the escape 
probability in treating emission lines. Dumont et al. use an "accel- 
erated lambda iteration" method that includes a detailed treatment 
of radiative transfer within the spectral lines themselves. They point 
out that their code has the drawback of being time-consuming, 
which precludes it from easily generating large grids of models for 
fitting data. In particular, they have applied their method to an iso- 
lated slab that is illuminated on one side. The total Thomson depth 
of the slab is tt = 4, and the illuminating flux corresponds to 
^ = 10"^ ergcms"^. As a test of our method, we have performed 
our own computation for this situation. We cannot expect complete 
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Figure 2. Spectra for the isolated slab illuminated on one side. The dashed 
curve is the illuminating spectrum, the solid curve is the "reflected" spec- 
trum emerging from the illuminated side, and the dotted curve is the "trans- 
mitted" spectrum emerging from the far side. The lower panel shows a 
close-up of a portion of the broader spectrum shown in the upper panel. 
The spectral resolution is 30. 



agreement because our calculations assume that helium is always 
fully ionized. Normally our interest is in reflection by a highly 
optically-thick medium. In order for all of the incident energy to 
escape back out the illuminated surface, the total radiative energy 
density increases with Thomson depth to such a degree that the as- 
sumption that helium is fully ionized should be satisfactory. When 
radiation is free to leave the far surface of a marginally thick slab, 
however, helium need not be fully ionized throughout. Dumont et 
al. (2003) found a large fraction of He^ ions at depths tt ^ 3. 
Our calculation of the isolated slab serves mainly as a test for the 
Lyman-Q lines of hydrogen-like ions, which are produced in the 
hotter layers nearer the illuminated surface. 

The biggest change required for this model is in the boundary 
condition at the far side of the slab. Diffuse radiation is allowed to 
emerge there, with the escaping flux calculated in the same way as 
at the illuminated surface (see Foster, Ross & Fabian 1986). Other 
changes must also be made to the method described previously. 
The gas density is reduced to nu = 10^^ cm^'^. The illuminating 
spectrum is a simple power law with photon index F = 2 extend- 
ing from 0.1 eV to 100 keV, and the energy resolution is fixed at 
E/AE = 30. The metal abundances are set to the values of Allen 
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(1973), which differ slightly from those of Morrison & McCam- 
mon (1983). 

Figure 1 shows the temperature structure and the ion fractions 
for oxygen and silicon that we calculate for the isolated slab. For 
tt < 2.5, our results are in good agreement with those of Du- 
mont et al. (2003). The differences at greater Thomson depth are 
presumably due to the lack of He"'' in our model. 

The reflected spectrum that we calculate for the isolated slab 
is shown in Figure 2. Dumont et al. (2003) concluded that use of 
the escape probability must result in a drastic overestimation of 
the strength of the Lyman-a lines of hydrogen-like ions. However, 
we find the Lyman-a lines of C VI (0.37 keV), O VIII (0.65 keV), 
Mg XII (1.5 keV), Si XIV (2.0 keV), and S XVI (2.6 keV) to be in 
excellent agreement with their results. In particular, the O VIII line 
has an equivalent width of 27 eV with respect to the total (incident 
plus reflected) continuum, compared to their value of 28 eV. For 
the Si XIV line, which is complicated by a steeper underlying con- 
tinuum and blending with the He-like line, we estimate an EW of 
14 eV compared to their value of 9.5 eV. In our calculations, Comp- 
ton scattering of resonance-line photons out of the narrow line core 
turns out to be much more important than reemission in the far line 
wings. That is, the first term of the numerator in equation (7) com- 
pletely dominates over the second term. 

The Fe Kq line is somewhat stronger in our calculation. Since 
it is dominated by the Fe XXV intercombination line (6.7 keV), 
however, this is probably due to different approximations in treating 
the atomic physics. Curiously, we find the Lyman-a line of N VII 
(0.50 keV) to be weaker than found by Dumont et al. (2003). Dif- 
ferences can be seen in other emission lines, but these are produced 
deeper within the slab, where He"'" could be a factor. Fig. 2 also 
shows the calculated spectrum emerging from the far side of the 
slab. Because of the lack of He"'", the EUV portion of the transmit- 
ted spectrum cuts off at a higher energy than found by Dumont et 
al. (2003). 

Overall, our method appears to give reasonably accurate re- 
sults, with computing times short enough to allow the production 
of a large grid of reflection models for fitting to observed X-ray 
spectra. 



3 RESULTS 

For an optically-thick medium, we have calculated model reflec- 
tion spectra covering a range of values for the ionization parameter 
(£, = 30, 100, 300, 1 000, 3 000, and 10000 erg cm s"^), the pho- 
ton index (P = 1.0 to 3.0 in steps of 0.2), and the iron abundance 
(^Pc = 0.1, 0.2, 0.5, 1.0, 2.0, 5.0, and 10.0) relative to its solar 
value. The result is a grid of 462 reflection models. 

Of course, the ionization parameter has a marked effect on 
emission and absorption features in the reflected spectrum. Figure 3 
shows reflected spectra for three different values of ^. In each case, 
the incident spectrum has F = 2.0, and iron has solar abundance 
(ylpc = 1). For ^ = lO** erg cm s^^, the surface layer is highly 
ionized, with Fe XXV-XXVII dominating for tt < 7. The only im- 
portant emission feature is the highly Compton-broadened iron Ka 
line peaking at 7.0 keV. For ^ — 10'^ erg cm s^^, the illuminated 
gas is not as highly ionized. Fe XXV dominates for tt < 1, and the 
lighter elements are no longer fully ionized below that. The strong 
iron Kq emission is dominated by the FeXXV intercombination 
line, and Compton broadening is still important. The spectral re- 
gion between 0.3 and 3 keV shows Kq lines of C VI, O VIII, Ne X, 
Mg XII, Si XIII-XIV, and S XV-XVI, as well as a few Fe Lq lines. 
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Figure 3. Reflected spectra for thi'ce values of the ionization parameter, 
with ^ = 10^ (bottom curve), 10'', and 10* erg cm s^^ (top curve). The 
incident spectrum has F = 2.0, and iron has solar abundance. 
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Figure 4. Reflected spectra for five values of the photon index, with F = 
1.2 (top curve), 1.6, 2.0, 2.4, and 2.8 (bottom curve). The incident spec- 
trum has t; = 10^ ergcms"^, and iron has solar abundance. Successive 
spectra have been offset by factors of five for clearer presentation. 



all atop a shallow absorption trough. The N VII Kq line, at an en- 
ergy (0.50 keV) just above the C VI K-edge, is very weak. When 
^ is reduced to 10^ erg cm s^^, there is very little ionization for 
'''T 5. The narrow Fe Kq line near 6.4 keV is due to fluores- 
cence, and the spectral region below 3 keV exhibits a myriad of 
emission features atop a deep absorption trough. 

In addition to affecting the overall slope of the reflected spec- 
trum, the photon index F also affects the emission and absorption 
features, since harder illuminating spectra have greater ionizing 
power. Figure 4 shows reflected spectra for five different values of 
F. In each case, the incident spectrum has ^ = 10^ erg cm s^^, and 
iron has solar abundance {Afc ~ 1). For F — 1.2 or 1.6, the illu- 
minated gas is more highly ionized than for F — 2.0. The iron Ka 
feature is a broader combination of Fe XXV and Fe XXVI lines, and 
the lines of the lighter elements are weaker. With F — 2.4 or 2.8, on 
the other hand, the illuminated gas is less highly ionized. Fe XVIII- 
XXII dominate at the outer surface instead of Fe XXV, so the Fe 
Kq line is much weaker, while the Fe Lq lines are stronger and 
more numerous. The Ka lines of the lighter elements are stronger 
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Figure 5. Reflected spectra for three values of the u'on abundance, with 
Ape = 0.2 (top curve), 1, and 5 (bottom curve) times solar abundance. 
The incident spectrum has § = 10'^ ergcms"^ and F = 2.0. Successive 
spectra have been offset by factors of five for clearer presentation. 

(as is the continuum absorption), with He-like exceeding H-like for 
Mg, Si and S. The oxygen lines are especially enhanced by a soft 
illuminating spectrum. 

The effects of altering the abundance of iron are shown in Fig- 
ure 5. Again, the illuminating radiation has ^ = 10'^ erg cm s^^ 
and r = 2.0. With iron reduced to 0.2 times its solar abun- 
dance, both the Ka emission line (6.7 keV) and the K-absorption 
feature (~ 10 keV) are reduced compared to the reflected spec- 
trum with solar abundance. Also, less Fe La emission is blended 
with the NeX Ka line (1.0 keV). With iron at five times solar 
abundance, the increased iron absorption lowers the continuum for 
1 < < 40 keV, while the iron Ka (6.7 keV), La (~ 1 keV), and 
2p 2s 0.1 keV) emission lines are enhanced. At the same 
time, the Ka lines of Mg, Si, and S are weakened, partly because 
these elements do not remain highly ionized to as great a depth. 

The grid of 462 models is being made available as an ionized- 
reflection model, called REFLION, for use in the XSPEC data-fitting 
routine (Arnaud 1996). 



4 DISCUSSION 

The treatment of Ross & Fabian (1993) has been extended to in- 
clude all commonly important ionization states and transitions. The 
resulting spectral grids should be useful for interpreting spectra 
from AGN, BHC and GRB. Major signatures of moderately ionized 
reflection are a Compton hump at about 30 keV, a strong ionized 
iron line at 6.7 keV, multiple emission lines on top of a soft contin- 
uum between 0.3 and 3 keV and a further EUV hump of emission 
peaking at about 60 eV. The iron abundance has a significant ef- 
fect on the shape of the reflected continuum between about 3 and 
30 keV. 

The spectra that we have calculated result entirely from re- 
flection (reprocessed illumination). A significant fraction of the in- 
cident X-ray energy emerges as softer radiation (as discussed re- 
cently by Madej and Rozaiiska 2004, for example). The accretion 
disc beneath the illuminated surface layer can be expected to pro- 
duce additional soft radiation (the "big blue bump") that can en- 
hance the EUV hump. Figure 6 shows what happens when the cal- 
culations include a soft blackbody spectrum with fcTbb = 0.010 
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Figure 6. Emergent spectra for pure reflection (bottom curve) and when a 
soft blackbody spectrum (fcT^b = 0.010 keV) enters the surface layer from 
below. Results are shown with the blackbody flux equal to the illuminating 
flux (middle curve) and ten times greater than the illuminating flux (top 
curve). The illuminating spectrum has ^ = 10'^ ergcms"^ and F = 2.0, 
and iron has solar abundance. Successive spectra have been offset by factors 
of three for clearer presentation. 



keV entering the illuminated layer from below. Models are shown 
with blackbody fluxes equal to the illuminating flux and ten times 
greater than the illuminating flux. In both cases, more EUV radi- 
ation emerges than for pure reflection, as expected, but the X-ray 
spectrum is essentially unchanged. Furthermore, accretion power 
liberated in an accretion disc corona (producing the illumination) 
should reduce the amount of power liberated locally within the disc, 
so the "big blue bump" may not be produced in the same regions of 
the disc as the X-ray reflection. For these reasons, we have chosen 
to calculate reflection spectra only. 

A signature of ionized reflection is the soft excess emission 
which occurs in the 0.2-2 keV band due to lines and bremsstrahlung 
from the hot surface layers. It gives a bump in the relativistically- 
blurred spectrum which, when folded though an XMM-Newton pn 
response matrix in order to simulate real data, is well-fitted by a 
blackbody of temperature 150 eV (residuals less than 10 per cent). 
This may be relevant to the 100-200 eV temperature component 
found in the X-ray spectra of many low-redshift PG quasars (Gier- 
linski & Done 2004; Porquet et al. 2004) which therefore may be 
due to ionized reflection. 

We have checked whether the many emission features in the 
reflection spectrum are likely to be readily observable. If there is lit- 
tle to blur the spectrum as a result of either instrumental resolution 
or shear within the ionized reflector, then the lines will be observ- 
able. However, when there is strong shear, as in the inner parts of an 
accretion disc, then the situtaion is much more difficult. We show 
in Figure 7 the effect of relativistic blurring in an accretion disc in- 
clined at 30 deg where the emission has an emissivity index of 3 and 
extends from 3 to 100 gravitational radii. Here ^ = 10'^ erg cm s^^, 
r — 2.0, and Apc — 1. The broad (helium-like) iron-K line is very 
clear, as is a relatively weak broad oxygen line. A weak feature of 
Fe-L and Ne with about the same equivalent width as O (~ 120 eV 
relative to the local reflection 'continuum' alone), together with Si 
and S are just discernible in this reflection-only spectrum. Figure 8 
demonstrates how extreme relativistic blurring and changes in iron 
abundance affect observations of iron-K features. A low iron abun- 
dance, perhaps due to SN II enrichment, produces a small edge in 
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Figure 7. Effect of relativistic blurring on the emergent spectrum. The 
model with ^ = 10^ ergcms"^, F = 2.0, and Ap^ = 1 (central line in 
Fig. 3, now dashed) has been blurred as if observed from a disc with emis- 
sivity index 3 extending from 3 to 100 gravitational radii around a maximal 
Kerr black hole, viewed at an inclination angle of 30 deg (solid line). 
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Figure 8. Effects of iron abundance on relativistic-blurred spectra. Models 
for 5 = 30 ergcms-i and T = 2.0 with Ape = 0.3 (solid) and = 3 
(dashed) have been blurred as if observed from a disc with emissivity index 
3 extending from 2 to 100 gravitational radii around a maximal Kerr black 
hole, viewed at an inclination angle of 30 deg. 

the spectrum, whereas a high iron abundance, perhaps due to SN la 
enrichment, produces an apparent large edge. 
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